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ABSTRACT 

We study the properties of voids in two different types of coupled scalar field theories. 
Due to the fifth force produced by the scalar field coupling, the matter particles feel 
stronger attraction amongst each other and cluster more quickly than they do in the 
standard ACDM model. Consequently voids in the coupled scalar field theories start to 
develop earlier and end up bigger, which is confirmed by our numerical simulations. We 
find that a significantly larger portion of the whole space is undcr-dcnscd in the coupled 
scalar field theories and there are more voids whose sizes exceed given thresholds. This 
is more prominent in early times because at later times the under-dense regions have 
already been evacuated in coupled scalar field theories and there is time for the ACDM 
model to catch up. The coupled scalar field theories also predict a sharper transition 
between voids and high density regions. All in all, the qualitative behaviour is different 
not only from the ACDM result, but also amongst specific coupled scalar field models, 
making voids a potential candidate to test alternative ideas about the cosmic structure 
formation. 
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1 INTRODUCTION 

One of the most active research areas in modern cosmology is 
about the theories involving cosmological scalar fields. As a 
potential candidate for dark energy (Copland et al 2006), for 
example, a canonic scalar field could have interesting proper- 
ties governed by its potential. Such quintessence (Wang et al 
2000) models are studied in depth in the literature, see- 
ing various potentials proposed, their properties investigated 
and their specific forms tried to be connected to develop- 
ments in high energy physics. Even richer phenomenology is 
achieved by considering variants of the simple quintessence 
model, such as giving the scalar field a non-canonical ki- 
netic term (Armendariz-Picon et al 2000) or coupling it to 
the matter fields (Amendola 2000) or even the spacetime 
curvature (Perrotta & Baccigalupi 1999), the latter case also 
covering certain modified gravity theories, such as the metric 
(Carroll et al 2005) and Palatini (Vollick 2003) f(R) grav- 
ity and Brans-Dicke theory (Brans & Dicke 1961). Often in 
these models, the behaviour of the scalar field is addition- 
ally controlled by a coupling function, and extra (fifth) forces 
are not uncommon. Local experiments then force the model 
builders to design specific mechanisms by virtue of which the 
(usually) severe experimental and observational constraints 
are not to be challenged, and remarkably these could be re- 
alised neatly by choosing appropriate scalar potential and/or 
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coupling function, as in the chameleon (Khoury & Weltman 
2004; Mota & Shaw 2007; Li & Barrow 2007; Hu & Sawicki 
2007) and environment-dependent dilaton (Brax et al 2010) 
models. 

An interesting property of these models is that either 
the mass of the scalar field or its coupling strength is sen- 
sitively dependent on the environments, in such a way that 
the fifth force is suppressed where the observations and mea- 
surements are made. Taking the chameleon model as an ex- 
ample: in high density regions the scalar field becomes very 
massive so that scalar field quanta could not propagate far 
and the fifth force gets suppressed. Exactly how massive the 
scalar field will be is determined by the local matter density 
as well as the steepness of the scalar potential and its deriva- 
tives, and the model could be so designed that the fifth force 
is suppressed in solar system but not on galactic and larger 
scales (Hu & Sawicki 2007), leaving the possibility that the 
structure formation process could be significantly affected. 

As one wants the scalar field mass to fluctuate strongly 
across the space, high degree of nonlinearity is inevitable. In 
chameleon model this is often reflected in the fact that the 
scalar field potential is very nonlinear, and in such circum- 
stances linear treatment obviously fails. A-body simulations 
are then the natural method to be used to study structure 
formation in these models, and in this paper we shall apply 
this very technique to study the void properties of them. 

Voids (van de Waygaert & Platen 2009) are an impor- 
tant ingredient of the standard picture of hierachical struc- 
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ture formation. Because the initial matter distribution in 
the Universe is inhomogeneous, as time passes by the over- 
dense regions will pull more matter towards them and under- 
dense regions get emptied to form voids. For the ACDM cos- 
mology, void phenomenon (Peebles 2001) are well studied 
(Tinker & Conroy 2009). In coupled scalar field models, the 
fifth force, if unsuppressed, will boost the clustering of mat- 
ter and therefore help to evacuate under-dense regions more 
quickly. A higher degree of emptiness in voids than what the 
concordance ACDM model predicts is thus an indicator of a 
possible fifth force (Keselman, Nusser & Peebles 2010). The 
voids are of even greater importance to chameleon models 
because, by naive expectation, the fifth force in these mod- 
els is suppressed in high density regions and shall not affect 
the galaxies clusters much, while in voids they are stronger 
and their effects more significant. It is therefore interesting 
to see what distinct features the voids have in these models 
compared with those in ACDM. As a preliminary work, we 
shall only consider dark matter voids (Colberg et al 2005) 
here, and leave the more technical work involving baryons 
to the future. 

In this work we will investigate two coupled scalar field 
models, one in which the scalar field is a chameleon and the 
other in which the scalar field is not. In practice, there is no 
sharp distinction between them. For the scalar potential de- 
scribed by equation (10) below, for example, the scalar field 
has no chameleon features if parameter a ~ O(0.1 — 1), but 
will become a chameleon while a C 1. Therefore, although 
for the non-chameleon models the fifth force is unsuppressed 
and thus could be well approximated as being proportional 
to gravity, as is commonly assumed by most A-body simu- 
lation works for coupled scalar field models to date (to name 
a few, Maccio et al (2004); Kesden & Kamionkowski (2006); 
Farrar & Rosen (2007); Hellwing & Juszkiewicz (2009); 
Keselman, Nusser & Peebles (2009); Baldi et al (2010)), for 
safety we shall solve the scalar field value as a function of 
spatial position explicitly and then differentiate to compute 
the fifth force (Li & Zhao 2009, 2010; Li & Barrow 2010a,b) 
(this technique has recently been applied to varying funda- 
mental constant (Li, Barrow & Mota 2010) theories as well). 

The arrangement of this paper is as follows: in § 2 we 
give a brief summary of the major ingredients of the coupled 
scalar field model. In § :? we explicitly write down the models 
we are simulating and present a detailed description of our 
void-finding algorithm. § 4 contains all our numerical results 
and finally § 5 is devoted to a summary and conclusions. 



2 THE BASIC EQUATIONS 

All the equations relevant for the simulations here are de- 
rived and discussed in (Li & Zhao 2009, 2010; Li & Barrow 
2010a) but for the present work to be self-contained we shall 
still list the minimum set of them which is necessary for us 
to understand the underlying physics. 

Instead of writing down the field equations directly as 
in some previous work, we start from a Lagrangian density 



K 



ip\/a<P 



+ V(ip) - C{ip)L OM + L s 



(1) 



in which R is the Ricci scalar, k = 8ttG with G being the 
gravitational constant, Ldm and Ls are respectively the 



Lagrangian densities for dark matter and standard model 
fields, ip is the scalar field and V(ip) its potential; the cou- 
pling function C(ip) characterises the coupling between ip 
and dark matter. Given the functional forms for V(ip) and 
C(ip) a coupled scalar field model is then fully specified. 

Varying the total action with respect to the metric g a b, 
we obtain the following expression for the total energy mo- 
mentum tensor in this model: 



T ab = Va^Vfc^ — gab 



V c V c ip-V{ip) 



+ C( ( p)T° M + T a s 



where T^ M and T^ b are the energy momentum tensors for 
(uncoupled) dark matter and standard model fields. The ex- 
istence of the scalar field and its coupling change the form 
of the energy momentum tensor, and thus could modify the 
cosmology from background expansion to structure forma- 
tion. 

Meanwhile, the coupling to scalar field produces a di- 
rect interaction (a.k.a. the fifth force) between dark matter 
particles, due to the exchange of scalar quanta. This is best 
illustrated by the geodesic equation for dark matter particles 



dt 2 



C(<p) 



V<p 



(3) 



where r is the position vector, t the (physical) time, $ the 
Newtonian potential and V is the spatial derivative. C v = 
dC /dip. The second term in the right hand side is the fifth 
force and only exists for coupled matter species (dark matter 
in our model). The fifth force also changes the clustering 
properties of the dark matter. Note that on very large scales 
ip could be considered as homogeneous and the fifth force 
vanishes. 

It has become obvious that in order to implement the 
above two equations numerically we need to solve both the 
time evolution and the spatial distribution of ip, and this 
could be done using the scalar field equation of motion 



V VaP H : h PDM : = 

dip dip 

or equivalently 

d(p 

where we have defined 
Veff(tp) = V(<p) + PdmC(vp). 



(4) 



(5) 



(6) 



The background evolution of ip can be solved easily once we 
know the current pdm, because pdm oc a~ 3 . We can then 
divide (p into two parts, ip = (p + Sip, where (p is the back- 
ground value and Sip is its (not necessarily small nor linear) 
perturbation, and subtract the background part of the scalar 
field equation of motion from the full equation to obtain the 
equation of motion for Sip. In the quasi-static limit in which 
we can neglect time derivatives of Sip as compared with its 
spatial derivatives (which turns out to be a good approxima- 
tion for our simulations, where the simulation box is much 
smaller than the observable Universe), we get 



^2 dCjip) 

V ip = PDM 

dip 



dC(tp) _ dV(ip) dV(<p) 
rz — Pdm + ■ 



dtp 



dip 



dp 



(7) 



where pdm is the background dark matter density. 

With the pdm made ready on some grid, we could then 
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solve Sip on that grid using a nonlinear Gauss-Seidel relax- 
ation method (in our simulations we have modified MLAPM, 
a public-available iV-body code using a self-adaptive refined 
grid so that high resolutions could be achieved in high den- 
sity regions). Because tp is also known, we then get the full 
solution of tp = (p+5ip. This then completes the computation 
of the source term for the Poission equation 

V 2 $ = 1 [C(ip) P DM - C(<f)pDM + Sp B - 25V(<p)] (8) 

where Sps = ps — Pb and 8V{ip) = V(tp) — V(p) are respec- 
tively the density perturbations of baryons and scalar field 
(note that we have neglected perturbations in the kinetic 
energy of the scalar field because it is always very small for 
our model). 

We could then solve equation (8) using a linear Gauss- 
Seidel relaxation method on the same grid to obtain $. With 
both $ and tp at hand, equation (3) could be used to com- 
pute the forces on the dark matter particles, and once we 
have the forces, we could do all the standard JV-body opera- 
tions such as momentum-kick, position-drift, time-stepping 
and so on. 

Equations (2 - 8) are all what we need to complete an 
JV-body simulation for coupled scalar field cosmology, and 
from then we could identify where the effects of the scalar- 
coupling come in: 

(i) The influence of the modified cosmic background ex- 
pansion rate mainly comes through the particle movements 
and time-stepping, i.e., equation (3). This is because in the 
simulations we shall use scale factor a as the time variable 
and d/dt = ad/ da. 

(ii) The varying mass effect could be seen directly from 
equation (8), which shows that the contribution of the dark 
matter density pDM to the source of the Poisson equation is 
multiplied by a factor C(ip) which differs from 1 in general. 
In our model the mass of dark matter particles is not really 
varying, but the net effect is just that. 

(iii) The fifth force appears explicitly in the right hand 
side of the geodesic equation (3), but only for coupled matter 
species (dark matter in our model). 

(iv) The velocity-dependent "frictional force" is a bit sub- 
tler. It hides behind the fact that equation (3), equation (7) 
are actually written in different gauges: equation (3) is the 
force for a dark matter particle and is given in that particle's 
rest frame, while equation (7) is written in the fundamental 
observer's frame. As a result, to use the Sip solved from equa- 
tion (7) in equation (3) we need to perform a frame trans- 
form VSip — > VS<p-\-a<pii in which x the comoving velocity of 
the said particle relative to the fundamental observer. This 
force is thus expressed as — jj-aipic, and obviously the faster 
a particle travels the stronger frictional force it feels. 

In our numerical simulation, we have included all these 
effects consistently (Li & Barrow 2010a). In particular, we 
have computed the fifth force explicitly, rather than simply 
assuming that it is always proportional to gravity: as shown 
in Li & Zhao (2009, 2010), such assumption could be fairly 
poor for certain models where the scalar field configuration 
is very inhomogeneous, although it is good enough for other 
models (Li & Barrow 2010a). 



3 SIMULATIONS AND VOID FINDING 
ALGORITHM 

3.1 The Models Studied 

In this work we consider the voids in the two different cou- 
pled scalar field models studied respectively by Li & Barrow 
(2010a) and Li & Zhao (2009, 2010); Zhao et al (2010). Both 
models have an exponential coupling between dark matter 
and the scalar field tp, 



C(tp) = exp^^/nip) 

and run-away potentials for ip: 

v(<p) = 



for the model of Li & Barrow (2010a), and 
V(<p) = 



exp 



(9) 



(10) 



(11) 



for the model of Li & Zhao (2009, 2010). In the above A is 
a parameter of mass dimension 4 and is of order the present 
density for dark energy [ip plays the role of dark energy in 
the models); its precise value is determined by the numer- 
ical code for the consistency in the background cosmology 
(Li & Barrow 2010a). 7 and a are dimensionless model pa- 
rameters controlling respectively the strength of the cou- 
pling and the steepness of the potentials. 

For the potential equation (10) we choose a = 0.1 and 
7 < so that the total effective potential V e ff(<p) is of run- 
away type. The scalar field then rolls quickly at early times 
((ptoday — Nearly ~ 0(Mp\) with A/pi the Planck mass), until 
slowing down to the slow-roll regime, when it behaves like a 
normal quintessence field. For the potential in equation (11) 
we choose a<l and 7 > 0, which ensures that V e // has a 
global minimum close to p = and d 2 V e f f(tp)/dtp 2 = m 2 at 
this minimum is very large in high density regions; then ip is 
trapped close to all through the cosmic history. These two 
cases are two extremes of the coupled scalar field field: in the 
former ip clusters very weakly just as in normal quintessence, 
and the fifth force is, to a good approximation, always pro- 
portional to gravity; the scalar field coupling also drastically 
modifies the background cosmology and structure formation 
at early times (z ~ 0(1 2 )). In the latter case ip is very inho- 
mogeneous and the fifth force is greatly suppressed in high 
density regions where tp acquires heavy mass, m 2 3> H 2 (H 
is the Hubble expansion rate), and thus the fifth force can- 
not propagate far. The suppression of the fifth force is even 
severer at early times, meaning that the structure formation 
is only influenced at late times (z less than a few); also, be- 
cause tp is trapped close to all the time, the background 
cosmology is forced to be indistinguishable from ACDM. In 
Table 1 we summarise the details for all the models we study 
here . 



For clearness we shall refer to the model with equation (11) as 
chameleon model, because here the scalar field mass m v depends 
sensitively on its environment and fluctuates strongly from point 
to point, and correspondingly the model with equation (10) as 
the non-chameleon model because here ra f is largely ignorant 
of the environment. The formal definition of a chameleon model 
(Khoury & Weltman 2004; Mota & Shaw 2007) is not relevant for 
our discussion in this work. 
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Table 1. A summary of the details of the models studied here. 
For all the runs: Q m = 0.257, n = 0.963, erg = 0.769, H = 71.9 
km/s/Mpc; size of simulation box is 64/i -1 Mpc and particle num- 
ber is 256 3 so that mass resolution is 1.04 X 10 9 /i — 1 Mg; domain 
grid has 128 cells in each side, and the finest refinement has 16384 
cells in each side, leading to a force resolution of ~ 12h~ 1 kpc. 
L, N, C stand for ACDM, non-chameleon and chameleon respec- 
tively. 



model 


v(v) 


C(ip) 


7 


a 




kA/3H% 


L 


constant 


1 


0.0 


0.0 




0.743 


Nl 


eqn (10) 


eqn (9) 


-0.05 


0.1 




0.717 


N2 


eqn (10) 


eqn (9) 


-0.10 


0.1 




0.781 


N3 


eqn (10) 


eqn (9) 


-0.15 


0.1 




0.838 


N4 


eqn (10) 


eqn (9) 


-0.20 


0.1 




0.891 


CI 


eqn ( ) 


eqn (9) 


0.5 


1.0 x 10 


-6 


0.743 


C2 


eqn ( ) 


eqn (9) 


0.5 


1.0 x 10 


—5 


0.743 


03 


eqn (11) 


eqn (9) 


1.0 


1.0 x 10 


-6 


0.743 


C4 


eqn (11) 


eqn (9) 


1.0 


1.0 x 10 


—5 


0.743 



3.2 Void Finding Algorithm 

Following Colberg et al (2005), our void finding algorithm 
consists of two steps: the identification of spherical proto- 
voids and mergers of proto-voids to form voids of arbitrary 
shape (all through this paper proto-voids and voids are dif- 
ferent things and are not to be confused). 

The proto-voids are spherical regions in which the aver- 
age of the density contrast 8 — p/p — lis below some prede- 
fined threshold S v . As shown by Colberg et al (2005), voids 
very clearly correspond to the troughs of the initial den- 
sity field, justifying the assumption that voids grow gravita- 
tionally from the initial negative overdensities. Assuming a 
spherical evolution model for the voids (Gunn & Gott 1972; 
Dubinski et al 1993), the growth of the voids can be studied 
analytically, and it is found that at the time of shell-crossing 
the overdensity inside the spherical proto-void reaches —0.8. 
Although this is the result for Einstein-de Sitter cosmology 
(Colberg et al 2005), we shall adopt it as a guidance and set 
S v — —0.8 in the coupled scalar field models as well. 

Our void finding algorithm is similar to that of 
Colberg et al (2005), but it differs from the latter in var- 
ious details, particularly the treatment of the mergers of 
proto-voids. To be clear and self-contained, here we briefly 
describe our algorithm in separate steps: 

(1) A regular 128 x 128 x 128 mesh is set up and 
the particle densities on this mesh are computed using the 
Triangular-shaped Cloud (TSC) scheme. This scheme ensures 
that the density interpolation is smoother than the usually 
used Cloud-in-Cloud scheme. 

(2) The local minima in the density field are located, 
and these are considered as the centres of the proto-voids 
(van de Waygaert & van Kampen 1993). Top-hat spherical 
windows with large enough radii so that the smoothed den- 
sity contrasts inside are greater than S v are then placed at 
these minima, and the radii are gradually decreased until the 
density contrast drops below S v . These minima and radii are 
then taken as the centres and sizes of the proto-voids respec- 
tively. One can also do this for all grid points on our mesh 
(i.e., set up a top-had window on each grid point, decrease 
the radii of the windows until the overdensity at a grid point 
falls below S v ) as in Keselman, Nusser & Peebles (2010), but 



this is more time-consuming and we have checked that the 
two methods lead to compatible results. 

(3) The above identified proto-voids are merged as ap- 
propriate to form the final voids of arbitrary shapes. It is 
well known that voids occupy the majority of the space, with 
islands of matter (dark matter halos and galaxies) intercon- 
nected by the narrow filaments which go through them. As 
a result, the merging criteria must be chosen carefully: for 
example, the dumbbell-shaped configurations are better to 
be avoided (Colberg et al 2005) to prevent the proto-voids 
from all being merged to form a single void as big as the sim- 
ulation box. To this end let's adopt a variant of the merging 
criterion proposed by Colberg et al (2005), which consists of 
the following steps: 

(i) Consider all the proto-voids which have not been as- 
sembled into any final voids. Find the biggest one, which is 
the primary progenitor of a merged void-to-be, and all the 
smaller ones which intersect it. 

(ii) If a smaller proto-void fully lies within the biggest 
one, then it is removed from the list. 

(iii) All the smaller proto-voids whose centres lie inside 
the biggest one are merged to the latter. 

(iv) All the smaller proto-voids whose centres lie outside 
the biggest one while still have significant intersection with 
the latter are merged to the latter. There are certain free- 
doms as for what is to be considered as significant, and here 
we will adopt the proposal of Colberg et al (2005): we divide 
the line segment connecting the centres of the smaller and 
biggest proto-voids into three sections, a section a which lies 
in both spheres and two sections b, c lying only in one of the 
two spheres respectively. The intersection is only considered 
to be significant if \a\ ^ max(jfcj, |c|). 

(v) All the smaller proto-voids which do not satisfy \a\ ^ 
max(|fo| , |c|) are not merged into the biggest one. In this case, 
the portion of volume shared by the smaller and the biggest 
proto-voids is assigned to the latter, and the volume of the 
smaller proto-void is decreased correspondingly. The smaller 
proto-voids are however not removed from the list but will 
be the building blocks for the voids-to-be considered later 
(as are the proto-voids which have no intersections with the 
currently biggest one). 

(vi) The total volume of the final merged void is that of 
the union of the biggest proto-void and all those smaller ones 
merged to (eaten by) it 2 . By our construction the shape of 
the final void could be arbitrary, but as Colberg et al (2005) 
we define an effective radius r e // such that ^§-r^ff equals the 
volume of it. The centre of the merged void, on the other 
hand, is taken to be the volume average of the centres of all 
merging blocks. 



2 The volume is computed as follows: we distribute a big number 
of points evenly in the whole volume so that the number density 
n is known, and then count the number N of points which lie 
inside the union. The volume of the union is then simply N/n. 
With our choices of the number density, we find that the numeri- 
cal results agree with analytical predictions (where the latter are 
available) typically better than 99.9999%. Wc are also careful to 
avoid assigning a same portion of space to more than one merged 
voids: for example, when computing the volume of a void-to-be, 
we always exclude the parts which lie in any of the previously 
identified voids. 
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(vii) The biggest proto-void and all the smaller ones eaten 
by it (but not the ones dealt with in step (v) above) are 
considered to be already assembled, and are excluded from 
the future runs of the merging process. Steps (i) - (vii) are 
then repeated until all the proto- voids are assembled. 

Note that using this algorithm we could naturally avoid the 
the dumbbell-shaped configurations for which two big proto- 
voids are connected by thin tunnels and are finally merged, 
because the two big proto- voids have no intersection by defi- 
nition, and are thus treated separately. The multiple merges, 
in contrast with the algorithm of Colberg et al (2005) which 
stops looking for further overlaps one it has found one, could 
produce bigger big voids, for which the biggest proto- voids 
typically eat tens of smaller ones, thereby expanding signif- 
icantly. 

In this work we only consider voids whose effective ra- 
dius r e ff 2.0h~ 1 Mpc, which is 4 times the size of the cells 
in the void-finder grid (Colberg et al 2005). As the majority 
of the empty space is populated by the very small voids, this 
will inevitably leave part of the void space undetected, as a 
price of maintaining accuracy. 



4 NUMERICAL RESULTS 

Having introduced the models and void finding algorithm in 
the above section, we now present and discuss the numerical 
results. 



4.1 Probability Distribution of Density Contrast 

As is shown in Li & Barrow (20I0a,b) and Li & Zhao (2009, 
2010), the existence of the scalar field and its coupling to 
dark matter efficiently enhance the clustering of matter, and 
we thus expect that an associated effect is the fast evacua- 
tion of the low density regions, resulting in a bigger portion 
of the space being inside voids (than in ACDM model). 

To give a first check of this expectation, we consider the 
probability distribution of the density contrast, which shows 
how much of the total space is in overdense/underdense re- 
gions. To obtain this, we place top-hat windows with radius 
2.0h~ 1 Mpc at each of our 128 3 regular grid points and com- 
pute the smoothed density inside them. We then count the 
number of grid points at which the density contrast S falls 
within [So, So + dS], and divide this by 128 3 to obtain the 
probability that S £ [So , So + dS] . 

The results for our 9 models are summarised in Fig. 1 
(see the figure caption for a detailed description), from which 
we could see clearly the trend that the coupling scalar field 
models (both N and C) predict more higher and lower den- 
sity regions than ACDM (L). However, there are some no- 
table differences between them, reflecting the fact that the 
N and C models are qualitatively different. 

For the C models, as mentioned above, the scalar field tp 
is trapped at the minimum of V e //, which is close to all 
through the cosmic history (as can be easily checked using 
the equation below). Using this fact it could be straightfor- 
ward to find that 
r- aA 

7PDM 

and thus (Li & Zhao 2009) 



2 

m v = 



d 2 V< 



dip 2 



cthiA 



(12) 



Although in reality the fifth force is quite complicated be- 
cause the force between two particles depends on the matter 
distribution in between them, the quantity m v nevertheless 
could be utilised to qualitatively understand the underlying 
physics. Basically, the greater m f is, the shorter the distance 
the fifth force could propagate, which means that a particle 
will feel the fifth forces exerted by less particles: in short, a 
heavy mass m v could suppress the fifth force. 

From equation (12), we see that m v increases as 7, Pdm 
increase or a decreases, and vice versa. Thus we expect the 
fifth force to be more severely suppressed in models CI, C3 
in which a is smaller, than in models C3, C4; within CI, C3 
it is more suppressed for C3 for which 7 is bigger, and for all 
models it is more suppressed at higher redshift, where pdm 
is bigger overall (Li & Zhao 2009, 2010). Furthermore, pdm 
is generally lower in the void regions, where the fifth force 
is expected to be less suppressed. 

In Fig. 1 we could see that at a = 0.5 (redshift 1, lower 
right panel) the probability distributions of the larger den- 
sity contrasts for the models C1-C4 are essentially the same 
as in model L, because fifth force is strongly suppressed there 
and then. On the other hand, for lower density contrasts in 
the low density regions, we do see the deviations from model 
L increasing for decreasing m v , as expected. As time passes, 
m v decreases overall and by a = 1.0 (redshift 0, upper right 
panel) a significant greater portion of space will be evacu- 
ated in the C models than in the L model, and meanwhile the 
enhancement in the clustering of matter due to the (attrac- 
tive) fifth force increases the chance the peaks in the initial 
density field develop into highly overdense regions (S 3> 1). 

For the N models, equation (12) does not apply because 
the scalar field always rolls down the effective potential V e /f. 
m 2 , does not fluctuate strongly in space, and the fifth force is 
never severely suppressed (Li & Barrow 2010a). As a result, 
the fifth force starts to affect the structure formation at a 
fairly early time. As is shown in Fig. 1, by a = 0.5 (redshift 
1, lower left panel) there have been significantly more over- 
dense and underdense regions than L model, mainly because 
the fifth force helps transfer more matter from low density 
to high density regions '. 

At a = 1.0, however, although the N models still pre- 
dict larger evacuated space than model L does (upper left 
panel), the volumes of space in the very overdense regions 
(S 3> 1) are less different between the N and L models, 
which is because as more and more matter is transferred to 
the high density regions, there is less and less remaining in 
the empty regions to be pumped: even though the aggre- 
gation of matter into high density regions starts earlier in 
the N models, it slows down eventually as matter in the low 
density regions is used up, and there turns out to be time for 
the L model to catch up somehow. These results for our cou- 
pled scalar field models (mainly the N models) is similar to 
that of the ReBEL model (Nusser, Gubser & Peebles 2005) 
as investigated by Keselman, Nusser & Peebles (2010) and 
Hellwing & Juszkiewicz (2009). 



3 Note that the probability is higher in the N models than in the 
L model for both l + i> 1 and 1 + S <C 1, as the total amount 
of matter is the same in both models 
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Figure 1. The probability distribution of the density contrast 8 (cquivalently 1 + 5). Upper Left Panel: Results at rcdshift for the 
non-chameleon models Nl (green dotted curve), N2 (blue dashed curve), N3 (purple dot-dashed curve) and N4 (pink dot-dot-dot-dashcd 
curve) in contrast to the ACDM (L) result (black solid curve). Upper Right Panel: The same but the green dotted, blue dashed, purple 
dot-dashed and pink dot-dot-dot-dashed curves now represent respectively the models CI, C2, C3 and C4. Lower Left Panel: The same 
as the upper left panel but for rcdshift 1. Lower Right Panel: The same as the upper right panel, but for rcdshift 1. 



4.2 Volume of Voids 

Having seen above that the fifth force in both the N and C 
models helps evacuate the low density regions, we now have 
a look at how the void properties are affected. 

The first interesting quantity is the void volume func- 
tion (similar to the halo mass function in the studies of the 
statistical properties of dark matter halos), which shows the 
number density of voids larger than a given volume V. Fig. 2 
displays our results for the N, C and L models. 

Understandably, the more strongly the matter particles 
cluster, the more effectively the low density regions are evac- 
uated and therefore the bigger the sizes of the voids tend to 
be. In the N models, not only does the fifth force, which is 
unsuppressed, start to take effect earlier, but also the uni- 
verse expands more slowly, leaving more time for particles 
to clump (Li & Barrow 2010a, b). Consequently, by a — 0.5 
(Fig. 2, lower left panel), we see large increases in both the 
void number density and void size in the N models as com- 
pared to the L model (for example, the number density could 
be twice as high). 

Again, as time goes on, low density regions are largely 



emptied and few particles remain there, then the growth of 
the voids will slow down. The earlier the evacuation starts, 
the earlier it will be completed and the growth of voids stops. 
As a result, when voids in the N models stop growing, those 
in the L model are still in the process. Finally, at a — 1.0, the 
void volume function of the L model more or less catches up 
those in the N models, as is evident in the upper left panel 
of Fig. 2. We also note that at a — 1.0 there are less small 
voids in the N4 model than in the L and other N models, 
which is likely because of the fact that small voids have been 
used up to merge to form bigger ones. 

For the C models, the suppress of the fifth force means 
that the clustering of matter and growth of voids are less af- 
fected by it 4 . This is easily seen in the a = 0.5 case (Fig. 2, 



4 One might think that, from equation (12), m v could be quite 
light in low density regions and thus the fifth force should be 
essentially unsuppressed there. The point, however, is that equa- 
tion (12) is only used to understand the physics intuitively, and 
in reality the fifth force is proportional to V<5</3 [cf. equation (3)] 
with Sip determined by the dynamical equation (7). As a result, 
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Figure 2. The void volume functions for the simulated models. Upper Left Panel: Results at redshift for the non-chameleon models Nl 
(green dotted curve), N2 (blue dashed curve), N3 (purple dot-dashed curve) and N4 (pink dot-dot-dot-dashed curve) in comparison to 
model L (black solid curve). Upper Right Panel: The same, but the green dotted, blue dashed, purple dot-dashed and pink dot-dot-dot- 
dashed curves now stand for model CI, C2, C3 and C4 respectively. Lower Left Panel: The same as the upper left panel but for redshift 
1. Lower Right Panel: The same as the upper right panel, but for redshift 1. In all the plots the horizontal axis is the void volume V, in 
unit of f/i _1 Mpc) , and the vertical axis is the number density of voids which are larger than V, in unit of lh~ 1 Mpc) 



lower right panel), which shows that the void volume func- 
tions for the C models do not deviate much from that for 
the L model (one might appreciate the effect of the suppress 
in the fifth force by considering that the ratio between the 
magnitudes of fifth force and gravity is 2^ if the former is 
not suppressed, and 7 ~ O(0. 1) for N models while 7 ~ 0(1) 
for C models). 

We could also have an examination of the void filling 
factor, defined as the fraction of total space that is filled by 
voids which are either bigger or smaller than V. Because our 
algorithm leaves the very small voids undetected, we choose 
to show the former, and the results are given in Fig. 3. It 
turns out that this plot shows more clearly the effects of 
the scalar coupling. As our first example, for the L model at 
a = 0.5 (Fig. 3, lower left panel), we notice that only 5% of 
the total space is filled by voids larger than 35/i -3 Mpc 3 , in 

the value of Sip in void regions generally depends on that outside 
the regions in the way that the solution to a differential equation 
depends on its boundary conditions. 



contrast to more than 11% and 18% for the models N3 and 
N4 respectively. At a — 1.0, as a result of void growth and 
mergers, the numbers for these three models are changed to 
25%, 35% and 45% respectively. In both cases, the scalar 
field coupling dramatically changes the total volume of void 
regions, and could potentially affect the properties of void 
galaxies (though the resolution of our simulations is limited 
and so we shall not touch this in this work). 

For the C models (Fig. 3, right panels), we obtain quali- 
tatively similar results, but the deviations from model L are 
obviously smaller due to the suppress of the fifth force. For 
example, at a = 0.5, the fraction of space occupied by voids 
larger than 35/i -3 Mpc 3 is respectively 5%, 6% and 7.5% for 
the models L, C3 and C4, while at a = 1.0 these numbers 
become 25%, 28% and 41%. 

4.3 Void Density Profiles 

The next quantity we are interested in is the void density 
profile, which characterises how matter is distributed within 
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Figure 3. The fraction of the total space that is occupied by voids larger than V as a function of V. Upper Left Panel: Results at 
redshift for the non-chameleon models Nl (green dotted curve), N2 (blue dashed curve), N3 (purple dot-dashed curve) and N4 (pink 
dot-dot-dot-dashed curve) in comparison to model L (black solid curve). Upper Right Panel: The same, but the green dotted, blue dashed, 
purple dot-dashed and pink dot-dot-dot-dashcd curves now stand for model CI, C2, C3 and C4 respectively. Lower Left Panel: The same 
as the upper left panel but for redshift 1. Lower Right Panel: The same as the upper right panel, but for redshift 1. In all the plots the 
horizontal axis is the void volume V, in unit of (h -1 Mpc) , and the vertical axis is the fraction of space occupied by voids larger than 
V. 



the empty regions. Like the growth rate of voids, this is also 
interesting and bears information about the physics driv- 
ing the structure formation. In Keselman, Nusser & Peebles 
(2010), for example, it is shown that the ReBEL model could 
produce different density profiles inside the void: the profile 
is steeper than the ACDM prediction, and is more so for 
smaller voids. In this subsection we would like to see what 
happens for our coupled scalar field models. 

As our voids are made from the spherical proto-voids 
and thus in principle could have arbitrary shape, it's difficult 
to give well-defined profiles for them. Instead, because we 
are only interested the steepness at the boundaries of the 
empty regions, here we choose to compute the profiles for the 
proto-voids, which are simply obtained by varying the radii 
of the top-hat smoothing windows located at the centres of 
the proto-voids and calculating the average densities inside 
them. 

As Keselman, Nusser & Peebles (2010), we'll consider 
two groups of the (proto-)voids, with radius ranges of 6 ^ 



To/ (h 1 Mpc) s£ 9 and 4 < r / (h 1 Mpc) s£ 6 respectively, 
where ro is the radius of the proto-void. 

For the group of larger voids, with radii in the range 6 ^ 
fo/ (/i -1 Mpc) ^ 9, the results are shown in Fig. 1. As in the 
ReBEL model, in both our N and C models the scalar field 
coupling makes the density profile deeper in the outer region 
of the proto-voids, and this effect is more prominent in early 
times (redshift 1.0, lower panels), which is not unexpected, 
because the fifth force helps to make the structure formation 
process start (and low density regions get evacuated) earlier, 
a fact which has also been confirmed by the observation that 
at early times the proto-voids have notably lower density in 
the very inner regions in the coupled scalar field models than 
in ACDM. Note that the chameleon effect does play a role 
here, by making the proto-void density profiles in models CI 
and C3 (for which a = 10 -6 and fifth force most severely 
suppressed) indistinguishable from that of the L model; at 
later times (redshift 0.0, Fig. 4, upper right panel), however, 
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Figure 4. The density profiles of the proto-voids whose radii fall in the range 6 ^ rg/ 1 Mpc) < 9. Upper Left Panel: Results 
at redshift for the non-chameleon models Nl (green dotted curve), N2 (blue dashed curve), N3 (purple dot-dashed curve) and N4 
(pink dot-dot-dot-dashed curve) in comparison to model L (black solid curve). Upper Right Panel: The same, but the green dotted, blue 
dashed, purple dot-dashed and pink dot-dot-dot-dashed curves now stand for model CI, C2, C3 and C4 respectively. Lower Left Panel: 
The same as the upper left panel but for redshift 1. Lower Right Panel: The same as the upper right panel, but for redshift 1. In all 
panels the horizontal axis is the distance from the proto-void centre, r, in unit of the radius of the proto-voids, rrj; the vertical axis is 
the dimcnsionlcss density contrast. 



the fifth force becomes less suppressed and the deviations 
from the L model gradually develop. 

Fig. 5 displays the same results, but for the group of 
smaller (proto-) voids. Here we can see the similar trend as in 
Fig. 4, namely the scalar coupling produces steeper increase 
in the void density profile around the void edge. The effect is 
not as strong as for the large halos, possibly because smaller 
voids form earlier and have been effectively evacuated even 
if the fifth force is not at play. 

These results suggest that larger voids at earlier times 
could better reveal the influences of a possible scalar field 
coupling than smaller ones. Also, the strong-chameleon and 
non-chameleon models can be distinguished because for the 
former the void profile is essentially the same as the ACDM 
prediction at earlier times but starts to deviate later, while 
for the latter the deviation starts quite early. 

Because of the limitation of resolution, we are unable 
to test the void profiles for very large and very small voids, 
which will be left for future work. Given the fact that large 



voids could be very different in size in different models, we 
expect that their density profiles could reveal more informa- 
tion about the physical models. 



5 DISCUSSION AND CONCLUSION 

Scalar-field-mediated long-range fifth forces have attracted 
much attention among cosmologists in recent years, and in 
most versions they come from a direct coupling between the 
matter species (usually dark matter only) and a cosmologi- 
cal scalar field. If they really exist, they might dramatically 
change the picture of cosmic structure formation, and allevi- 
ate or even solve some problems in the concordance ACDM 
model. On the other hand, the range and magnitude of the 
fifth force are often model dependent: in one extreme, rep- 
resented by our N models, the effective potential V e j / (<p) is 
fairly flat so that the mass m v is light and almost the same 
everywhere; the fifth force then has a fixed ratio to grav- 
ity, which is equivalent to a rescaling of the gravitational 
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Figure 5. The same as Fig. 4, but for proto- voids whose radii fall into the interval 4 ^ ro/ (h 1 Mpcj < 6. 



constant. In the opposite extreme, like our C models, V e ff 
and thus the scalar field mass m v depends sensitively to the 
matter density in the way that m v could be very heavy in 
high density regions, where the fifth force is severely sup- 
pressed and thus negligible, but is quite light in low density 
regions, where the fifth force has a fixed ratio to gravity as in 
the N models. Furthermore, as emphasized by Li & Barrow 
(2010b), the fifth force is often not the only impact the cou- 
pled scalar field could have on cosmology, nor even is it usu- 
ally the most important one. In the N models, for example, 
the modification of the cosmic background expansion rate by 
the scalar coupling could be more influential in the course 
of structure formation (Li & Barrow 2010b). 

The complexities indicate that the model-independent 
studies of the coupled scalar field might fail to account for 
the various effects due to the scalar field coupling (see § 2 for 
a description) appropriately. So in this paper, we have stud- 
ied the formations and properties of voids in the L, N and C 
models in parallel and compare their predictions. Voids are 
the largest objects in the Universe which are produced dur- 
ing the course of structure formation, and fill the vast ma- 
jority of the space. The importance of their properties in un- 
derstanding the underlying cosmological scenario and global 
cosmological parameters has been emphasised by many au- 
thors. Recently, it has been claimed that a long-range fifth 



force could evacuate the space more efficiently and thus pro- 
duce more voids than ACDM (Keselman, Nusser & Peebles 
2010; Hellwing & Juszkiewicz 2009). Those studies concen- 
trate on the ReBEL model, where only a Yukawa-type fifth 
force is considered; because of the reason mentioned above, 
here we take into account all the main effects due to a scalar 
field coupling (which is arguably the most natural cause of 
a fifth force), consider two qualitatively different types of 
models and make more detailed analysis of their effects on 
void properties by revising the void-finding algorithm. 

Our Fig. 1 shows that in the coupled scalar field mod- 
els matter is more concentrated in some regions, leaving the 
remaining of the space more evacuated, in line with the ex- 
pectations. Here we note that the N and C models behave 
differently: for the former, the fifth force is never suppressed, 
and the migration of matter from low density regions to high 
density regions starts earlier thanks to it; for the latter, the 
chameleon effect suppresses the fifth force at the early times, 
the effect of which in boosting the clustering of matter only 
becomes significant recently (after redshift 1). In both mod- 
els, larger portion of space is under-dense today than in the 
L model. 

We apply our void-finding algorithm to the N and C 
models, and find that both models predict a bigger number 
of larger voids than the L model does (Fig. 2). Once again, 
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the N and C models behave quite differently, in particular at 
early times: by redshift 1.0 the N4 model produces several 
times more voids than the L model (and also the biggest 
voids are several times bigger), while the C models are only 
slightly different from the L model, though the fifth force in 
them, if unsuppressed, is much stronger than in the N mod- 
els. The result seems to be contradictory to the expectation 
that in the void regions the fifth force gets less suppressed 
and therefore we should have seen greater difference from 
the L model. The reason is as follows: firstly, the scalar field 
equation of motion is dynamical and the solution in the void 
regions depends on the overall environment in the simula- 
tion box, so it is untrue that the fifth force in void regions is 
unsuppressed; secondly, the formation rate of voids is more 
dependent on how fast high density regions could pull mat- 
ter out of them, but as the fifth force is suppressed this pull 
is not much stronger than in the L model. 

It is worth noting that the difference between the void 
volume functions in N and L models is bigger at earlier times 
(Fig. 2), because at later times most potential void regions 
have already been developed: since there are not many more 
new voids yet to be produced in the N models, the L model 
gradually catches up. For the C models, the trend is quite 
opposite, and more voids are produced recently than in the 
L model, because finally the fifth force is freed and starts to 
take effect. 

We also show in Fig. •'! the fraction of space which is 
filled by voids exceeding a certain threshold in volume. Here 
the qualitative features could be explained by the same ar- 
gument used for Fig. 2, and what is more impressive is the 
quantitative result it illustrates. For example, at redshift 1.0 
the N4 model, which is the most extreme in the N models, 
predicts almost 4 times as much space filled by voids larger 
than 35 /i _3 Mpc 3 as does the L model. Even at present the 
number is almost 2 (the same also applies for the C models) . 
Voids prove to be a promising tool to constrain the scalar 
field couplings. 

Finally, we have studied the density profiles of the voids 
(Figs. 1, 5). We find that in general the voids in the coupled 
scalar field models are featured by a sharper transition from 
low density to high density around their edges, similar as the 
result for the ReBEL model (Keselman, Nusser & Peebles 
2010). At earlier times the large voids in coupled scalar field 
models also have lower overdensities in the inner part due 
to the more effective evacuation of the region. 

Due to the limitation of the simulation box and resolu- 
tion, we have not studied other interesting void properties 
such as the halos in voids. And because we have only dark 
matter in the simulations, we have not touched the forma- 
tions and properties of the void galaxies. These shall be left 
to future works. The existing results, however, already indi- 
cate that the void properties could be largely influenced by 
a coupled scalar field, being it chameleon or not, and there- 
fore voids could provide a useful tool to study and constrain 
such alternative scenarios for structure formation. 
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